Coherent Inverse Photoemission Spectrum for Gutzwiller Projected Superconductors 
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Rigorous relations for Gutzwiller projected BCS states are derived. The obtained results do not 
depend on the details of model systems, but solely on the wave functions. Based on the derived 
relations, physical consequences are discussed for strongly correlated superconducting states such 
as high-Tc cuprate superconductors. 
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Right after the discovery of high-Tc cuprate super- 
conductors , Anderson has proposed a Gutzwiller pro- 
jected BCS wave function — a quantum many-body state 
incorporating strong on-site Coulomb repulsion — to de- 
scribe the superconducting state 0]- Since then there 
have been extensive studies in understanding the nature 
of this state and its variants 0: 0|- In addition, sev- 
eral reports have shown that this projected BCS wave 
function is indeed a good variational ansatz state to de- 
scribe the ground state of t-J hke models 0, S H 0' 
which are believed to capture the low energy physics of 
the cuprates 0. Although these projected BCS states 
were proposed more than 15 years ago, very recently 
they have acquired a revived interest 0, 0, IT3| . 
This is probably because recent expensive numerical cal- 
culations based on the Gutzwiller projected variational 
ansatz clearly indicate that many aspects of the physics 
of high-Tc cuprate superconductors can be understood 
within this framework p^ . 

In this short communication some rigorous relations 
are derived for the Gutzwiller projected BCS states. It 
is shown that, as a consequence of the derived relations, 
the one-particle added excitation spectrum tends to be 
more coherent than the one-particle removed excitation 
spectrum docs. It is further shown numerically that this 
trend is still observed approximately for more involved 
Gutzwiller projected BCS states. Possible experimental 
implications of the present results are also discussed. 

Here our general system consists of a single orbital per 
unit cell on the two-dimensional (2D) square lattice with 
L sites 01 . The creation and annihilation operators of 
spin a{=1, i) particle at site i are denoted by C;^ and Cia, 
respectively. A Gutzwiller projected BCS state with N 
particles is described by 



If^^') -P^PgIBCS), 



(1) 



where P^r is the projection operator onto the fixed num- 
ber N of particles, Pq =Y[ii^~ "-iT^ii) ^^"3 Gutzwiller 
projection operator to restrict the Hilbert space with 
no double occupancy on each site, and fii(j — Cj^Cjcr. 
|BCS) — Hk a 7kcr|0) is the ground state of the BCS mean 
field Hamiltonian where 



are the standard Bogoliubov quasi-particle operators, 
Cka — X]i'2^'''''cicr/\/T, |0) is the vacuum of particles, 
and the singlet pairing is assumed T^. The nature of 
this state has been extensively studie d sp ecially in the 
context of high-Tc cuprates |1 H ItL 13. ITo| 

A one-particle added state with spin a and momentum 
k is similarly defined by using 7^^: 



(3) 



This state was first proposed by Zhang, et al , followed 
by several others [lol u3 . [Til l . Hereafter the normalized 
wave functions for the N- and {N + l)-particle states are 

denoted by IV'o^'') and IV'k^'^^^)' respectively. 

First it is useful to show that the following operator 
relation between Cko- and cj^.^ holds exactly: 



kcr 



Pg 



Ckcr-PcC, 



kcr 



G- 



(4) 



Here No- = X^i ^L^i'^ ^^'^ ^ stands for the opposite spin 
of a. This is easily proved by using Pg<^s.„Pg — -PgC;^- 

Using the above equation it is readily shown 

that the momentum distribution function no-(k) — 



(i/'o^''|cj^^Ckcr|'0o™'') calculated for the state l^'g"'') is re- 



lated to the state l^" 



kCT 



through 



(Ar+l)l^(7V+l) 



Wk 



k(T 



ko- 



(5) 



where TSS„Pq — PcT^a is used. 

The quasi-particle weight for the one-particle added 
excitation is defined by 



(+) 



(6) 



Now we shall show that there exists a simple and ex- 
act relation between z'^^ and rto-(k). To this end, it is 
important to notice that since 



7kT 
7lk, 



"k -"^k 

Mk 



CkT 



(2) 



|5LI*r)-"k( 



kV< 



+l)l^(w+l) 



k(T 



(7) 
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is simplified as 



(+) 



(8) 



From Eqs. (|SJ and (O, we finally arrive at the desired 
relation, 



(9) 



It should be emphasized that to derive the above equa- 
tion we have not made either any approximations or any 
assumptions except for the form of the wave functions for 
the N- and {N + l)-particle states given by Eqs. iQJ and 
©, respectively. 

The equation is also simply verified numerically on 
small clusters using a Monte Carlo technique. Typical 
results are presented in Fig. ^ As seen in Fig. ^ and 
for all other cases studied, Eq. @ is satisfied within the 
statistical errors. 
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FIG. 1: no-(k) (circles) and Z^^^ (squares) calculated using 
a Monte Carlo technique for L = 16 x 16 and N^- = — 
115 |l^ . The sum of the two quantities {na(k) + Zl^^) is also 
plotted by crosses, which are 1 — N^/L = 0.551 within the 
statistical error bars (smaller than the size of the symbols). 



In order to discuss a physical consequence of Eq. ^ on 
the one-particle excitation spectrum, let us first derive a 
simple same rule. The one-particle excitation spectra for 
removing one particle [A^^^(k, w)] and for adding one 
particle [^J^^^^(k, cj)] are defined respectively by 



Ar^(k,c.) = --Im/cL 



1 



uj + {H - Eo) + iO+ 



st 



TT \ U! - (H - Eq) + lO+ / 



where (O) is the expectation value of O for the exact 
ground state of the system described by Hamiltonian H 
with its eigenvalue Eq. It is generally proved that for 
any model systems defined within the restricted Hilbert 
space with no double occupancy on any sites by particles, 
such as the t-J model, the 0-th moment of the spectral 



function satisfies the following sum rules: 

/oc 
-OO 



(10) 
(11) 



where (k, = Al^^{k,u) -f A^^^^ [k, uj) . The lat- 

ter equation is easily proved by using Eq. Eq. ifTUI) 
is a rather standard sum rule, while Eq. (|ll|l is due to 
the reduction of the Hilbert space by Pq and it is indeed 
satisfied for, e.g., the t-J model [l8l |. 

Let us now discuss what physical consequences would 
be expected. First we assume that there exists a sys- 
tem for which the ground state and the low-lying excited 
states are approximately described by the wave functions, 
Eqs. and (j2Jl, introduced above. Then it follows im- 
mediately from Eqs. (P)- (fTT|l that 



A^/^==(k,w)dw = Z, 



(+) 

kcr ■ 



(12) 



This relation implies that the one-particle added excita- 
tion spectrum is all coherent since only one state con- 
tributes to A}^^^{k,Ljj). It should be noted here that 
while several studies have recently reached the similar 
conclusions |l2lll3j |. the argument presented here is more 
rigorous and transparent. 

Next we shall discuss to what extend Eq. holds and 
therefore Eq. (|12|l remains approximately true for more 
involved wave functions. A natural and important ex- 
tension of the simplest Gutzwiller projected BCS states 
described by Eqs. ^ and ||2Jl can be achieved by includ- 
ing charge Jastrow factors Jc, *.e.. 



1$, 



FatFgJcIBCS), (13) 
\'^\i7") = Pn+iPgJc^JBCS), (14) 



'0 / 



for the ground state and the one-particle excited states, 
respectively. Here 



Jc = exp 



(15) 



ni = nil -|- fiii , and the sum runs over all the independent 
pairs of sites i and j. The importance of Jc has been 
already reported for various lattice models [3, A 
typical example of Vij is presented in Fig.|2(a) where all 
the independent are optimized for the 2D t-t'-J model 
with J/t = 0.3, t'/t = -0.2, and ^ ^ 115 on 
L = 16 X 16 [IQ]. Certainly the inclusion of Jc improves, 
e.g., the variational energy. Besides such quantitative 
changes, Jc can also make a qualitative difference. One 
of these examples is shown in Fig.|2](b), where the charge 
structure factor A^(q) = exp(— iq- l)(nini+i) for small 
wave numbers are calculated using the wave functions 
with and without Jc. As seen in Fig. [21(b), A^(q) 
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as |q| for |<&g^^) as expected, whereas iV(q) finite 

as |q| ^ for l^-^,^') This is because |BCS) does 

not conserve the number of particles which is instead a 
conserved quantity for the t-t'-J model [2ll[2^ . 



So far we have only considered the one-particle added 
excitations. Let us briefly discuss the one-particle re- 
moved excitations. A one-particle removed state is anal- 
ogously constructed by 
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I* 



(JV-l)v 



Pn-iPgiL\^CS). 



(16) 



Ahhough 1*^,^"^' 



and 1* 



'k(T / ^'"^ very like form, 
the similar conclusion about the coherence of the one- 
particle excitations can not be drawn for the one-particle 
removed excitations. A simple reason for this is the fol- 
lowing 12, 13]: for the one-particle added excitations, 



FIG. 2: (a) Charge Jastrow factor v{r) — uij as a function of 
distances r = |i — j[. These quantities are optimized in such a 
way that the variational energy of I^'q^'') is minimized for the 
2D t-t'-J model with J/t = 0.3, t'/t = -0.2, and iV| = iV| = 
115 on L = 16 X 16 (b) Charge structure factor iV(q) 

calculated using l*^,'^') (circles) and (crosses) for the 

2D t-t'-J model with the same model parameters as in (a). 
The variational parameters are optimized for both states. The 
statistical error bars are smaller than the size of the symbols. 
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FIG. 3: no-(k) (circles) and (squares) calculated using 

the wave functions |<I>[,^') and with Jc- The model 

parameters used are the same as in Fig. |5| The sum of the 
two quantities {ncr{ii) + Z^^) and 1 — N^/L = 0.551 are also 
plotted by crosses and thick dashed lines, respectively. The 
statistical error bars are smaller than the size of the symbols. 



It is now interesting to examine if the exact relation 

Eq. © proved for the states |*[,^^) and 1*^,^"^^^), and 
thus Eq. (|12() , can still hold for these more involved wave 

functions I'&o^'') and I'&k^^^''). The numerical results on 
small clusters are presented in Fig 13 for the 2D t-t'-J 
model with the same model parameters as in Fig. |21 As 
seen in Fig 13 surprisingly Eq. Q remains satisfied, at 
least approximately. Thus it can be still argued that, 
because of the sum rule for the one-particle excitation 
spectrum [Eqs. ((Tn|l and ifTT)) ]. Eq. ((T^ is approximately 
satisfied, and therefore most of the one-particle added 
excitation spectrum consists of a coherent part. Further 
numerical calculations have been carried out for the 2D 
t-t'-J model with different model parameters, and it was 
found that Eq. ® is stiU satisfied within 10-15% 



|BCS) (X |vl/(^+i)) 



{N)\ 



consists of only one state, while for 



the one-particle removed excitations ^c'Ckcrl^'o 



(N) 



, which is not described by l^'^L^^^'') 



alone. It is 



also checked numerically on small clusters that the quasi- 
particle weight for the one-particle removed excitation, 



/ 



is substantially different 



from n(j(k)- 

Finally we shall discuss experimental implications of 
the present results. A most relevant experiment is angle- 
resolved inverse photoemission spectroscopy on the su- 
perconducting state for the hole-doped cuprates. If we as- 
sume that the Gutzwiller projected BCS states discussed 
here are faithful description for the superconducting state 
in the cuprates, it is expected that the inverse photoemis- 
sion spectroscopy spectrum has more coherent character- 
istics than the direct photoemission spectroscopy spec- 
trum does [231 . The similar trend is also expected in 
the superconducting state for the electron- doped cuprates 
except now that the direct photoemission spectroscopy 
spectrum has more coherent characteristics. This is be- 
cause the t-J like models can also describe the electron- 
doped cuprates only after the particle-hole transforma- 
tion: Cko- — > ^k+Qcr [Q = 113 J and the same argu- 
ment presented here is still true for h'^^. 

To summarize, we have derived some rigorous rela- 
tions for the Gutzwiller projected BCS states. Using a 
sum rule for the one-particle excitation spectrum, it was 
shown that the one-particle added excitation spectrum 
tends to be more coherent than the one-particle removed 
excitation does. Possible experimental implications were 
also discussed. Finally, it should be noted that all the 
results presented here are based on the Gutzwiller pro- 
jected BCS states studied, and a question of whether 
these states can represent the exact eigenstates of some 
particular model Hamiltonian is beyond the scope of the 
present study. 
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